Phase transitions in Ising model on a Euclidean network 
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A one dimensional network on which there are long range bonds at lattice distances I > 1 with 
the probability P(l) oc l~ 6 has been taken under consideration. We investigate the critical behavior 
of the Ising model on such a network where spins interact with these extra neighbours apart from 
their nearest neighbours for < S < 2. It is observed that there is a finite temperature phase 
transition in the entire range. For < 5 < 1, finite size scaling behaviour of various quantities are 
consistent with mean field exponents while for I < S < 2, the exponents depend on 8. The results 
are discussed in the context of earlier observations on the topology of the underlying network. 
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I. INTRODUCTION 

During the last few years there has been a lot of ac- 
tivity in the study of networks once it was realised that 
networks of diverse nature exhibit some common features 
in their underlying structure. The Watts-Strogatz (WS) 
model [J was proposed to simulate the small-world fea- 
ture of real networks (i.e., the property that the networks 
have a small diameter as well as a large clustering coeffi- 
cient). In this model, the nodes are placed on a ring and 
each node has connection to a finite number of nearest 
neighbours initially. The nearest neighbour links are then 
rewired with a rewiring probability p r to form random 
long range links. This model displays phase transitions 
from a regular to a small world to a random network 
by varying the single parameter p r representing disor- 
der. The regular to small world network transition was 
shown to take place at p r — > 0. In a slight variation of 
the WS model, the addition type network was consid- 
ered in which random long range bonds are added with a 
probability p, keeping the nearest neighbour links undis- 
turbed. Phase transition to a small world network was 
again observed for p — * 0, ■ 

Phase transitions in networks can also be driven by fac- 
tors other than the above kind of disorder. In many real 
world networks, the linking probability is dictated by fac- 
tors like Euclidean distances separating them, aging etc. 
It is possible to achieve phase transitions in the network 
by tuning the parameters governing such factors. These 
networks are indeed not just theoretical concepts. Many 
real world networks like the Internet at the router level, 
transport networks, power grid network and even col- 
laboration networks are indeed described on a Euclidean 
space in which the geographical locations of the nodes 
play an important role in the construction of the net- 
work US (10. Similarly, aging factor is important in 
social networks and citation network 0, 0, where it 



is found that linking with older nodes are generally less 
probable. 

In the theoretical modelling of Euclidean networks E3 , 
it is usually assumed that two nodes separated by a Eu- 
clidean distance I are connected with the probability P(l) 
which follows a power law variation, i.e., 
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p(i) ~ r 
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In one dimension, where the nodes are placed on a ring, 
the typical networks generated for different values of 
5 and the corresponding adjacency matrix are shown 
graphically in Fig. ^ Growing networks on Euclidean 
spaces where the linking probability is modified with such 
a probability factor have also been considered recently 
jl 1| . In a real situation, however, P(l) may not have 
such a well-defined behaviour [1 Q • 

A variation of P(l) in the above form was first stud- 
ied by Kleinberg [l2| on a two dimensional plane with 
the aim to find out how navigation in the network 
depends on the paramater 5. Later, a number of 
other properties of networks (both static and dynamic) 
with such a probabilistic attachment have been studied 

niEiEiEiEaEa 

In the present study, we have reconsidered a static one 
dimensional Euclidean network (OEN) with connection 
probability at distance I given by eq We are inter- 
ested in the phase transitions which can be achieved by 
varying the parameter 5. Although this issue has been 
addressed earlier 0,0,0,0], there remain some ques- 
tions about the exact nature of phase transitions in the 
region < 5 < 2. We have tried to resolve the problems 
by taking a different approach here. In the earlier stud- 
ies, the topological features like the diameter, shortest 
paths and the clustering coeffcient of the network had 
been studied. Here we have considered Ising spins on 
the nodes of such a network. The critical behaviour of 
this system is expected to reflect the nature of the net- 
work at different values of S indirectly. That the Ising 
model undergoes phase transitions on the addition type 
WS network is an established fact 0,0, |2(j- It has also 
been shown that in the small world pha se the addition 
type WS has a mean field nature |2J, |23 ■ 
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In section m we have reviewed the available results of 
the Euclidean networks and also justified the necessity 
of the present investigations. In section IIIII studies of 
Ising models on networks have been briefly reviewed. Our 
results of the study of the Ising model on a Euclidean 
network is given in section IIVI Finally, in section 
summary and conclusions are given. 



II. NETWORK ON EUCLIDEAN SPACE 

The one dimensional Euclidean network (OEN), with 
the connection probability given by eq l[T|l. is equivalent 
to the addition type WS model at 5 = when addition of 
long range links can take place at any length scale. The 
choice of the number of long range neighbours in this 
OEN should be such that it allows a phase transition at 
(5 = 0. In the addition type WS model, this transition 
for a system of N nodes is achieved for p = 1/N (im- 
plying p — > in the thermodynamics limit). Connection 
probability p = 1/N implies a total number of N long 
range edges in the system or the existence of one long 
range bond per node on an average. Therefore in order 
to obtain a phase transition in the OEN, it is sufficient 
to keep N long range bonds. For large values of 5 the 
nodes in the network have short range connections only 
and therefore it behaves as a regular network 0,0,^3- 
Thus there should exist at least one phase transition in 
this network. 

In all, there could be four kinds of behaviour of this 
network: (a) Regular: when the network is like a short 
ranged one dimensional system (b) finite dimensional: 
when the network behaves as a system with effective di- 
mensionality greater than one but still short ranged (c) 
Small world and (d) Random. A few studies in which 
these phase transitions have been investigated have not 
been able to give a concrete or unique answer. The net- 
work behaves as a regular network above 6 — 2; this is 
obtained in all the earlier studies. The interesting region 
is < S < 2 as here the network no longer behaves like 
a one-dimensional system. In |13). the behaviour of the 
OEN was studied by releasing a random walker on the 
network. The results indicated that the walker has the 
same behaviour as on a small world network for all 6 < 2. 

To detect whether a network has small world be- 
haviour, one can calculate the shortest distances separat- 
ing two nodes, take the average and analyse its behaviour 
with the system size. The largest of these shortest paths 
is called the diameter of the network. The average short- 
est path and diameter of a network are expected to have 
the same scaling behaviour. Sen and Chakrabarti 0] 
studied the diameters of the OEN while the shortest 
paths were calculated by Moukarzel and de Menezes 0] 
giving contradictory results. In jlj], it was found that 
the diameter behaved as In N on chains of size N for all 
5 < 2 and hence it was concluded that small world be- 
haviour occurs for < <5 < 2 while in it was argued 
that small world behaviour occurs only for S < 1. In the 




site i 

FIG. 1: Structures of the networks (left panel) and the corre- 
sponding adjacency diagrams (right panel) for <5 = 0, f .0, 2.0 
(top to bottom) in the model for a system of size N = fOO. 



region 1 < S < 2, according to 0, the shortest distances 
scale as N with the value of the shortest path exponent 
6 (0 < 6 < 1/2) depending on <5. The system sizes con- 
sidered in were much larger compared to those in 0] 
and the discrepancy of the results could be accounted for 
by this fact. It must also be mentioned here the above 
results are for pN — 1 , results for other values of p have 
also been considered in these studies. 

Random and small world phases can be distinguished 
by the clustering coefficient which remains finite in the 
latter. The estimate of the clustering coefficients of the 
present network |l6j | showed that below S = 1 it vanishes. 
This was considered to be a signature of the network be- 
ing random for S < 1 and based on the findings of [l4| . 
it was concluded that the region < <5 < 2 was equally 
divided into a random (for 6 < I) and a small world (for 
5 > I) phase. But this is by no means foolproof as in this 
particular network the clustering coefficient is bound to 
be small close to 5 — as the number of nearest neigh- 
bour is small when S is small. (In the original WS model, 
the clustering coefficient remained high for finite values 
of p r ^ 1 only when the number of nearest neighbours 
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FIG. 2: Behavior of shortest path S with system size N. 



was at least four to begin with.) Thus it may not indi- 
cate a random to small world phase transition at 5 = 1 as 
conjectured in [iHf. it could, in principle, correspond to 
the small world to a finite dimensional network as found 
in [l^. That there could be two transitions, one at 5 = 1 
and the second at 6 — 2 was also supported by simple 
scaling arguments for the average bond lengths [lfij . 

It has been mentioned in the beginning of this section 
that the the total number of random long range bonds 
per site has been kept equal to 1 on an average such that 
p = 1/N. This is also consistent with the fact that such a 
choice will keep the probability normalised. However, in 
an analytical approach there will be a problem at S = 1 
where such a normalisation is not possible. This issue 
has been discussed in [l3j and does not cause a problem 
in numerical calculations. 

Navigation or searching on small world networks are 
known to give rise to shortest paths which do not behave 
as ln(iV) but rather have sublinear variations with N 
except for special points: this was first detected in 
for two dimensional lattices and later confirmed for the 
WS model |23 . In the present one dimensional Euclidean 
model this was again demonstrated in [l7| . 

The more recent results indicate that there is 

indeed a phase transition at S = 1 and that the network 
behaves as a small world only for S < 1. In the present 
work we want to confirm this using a completely different 
approach. If the network's behaviour is that of a small 
world it should be reflected in the critical exponents of 
the Ising model which would assume mean field values. 
On the other hand, if it behaves like a finite dimensional 
system with effective dimension greater than one, a phase 
transition will be observed with critical exponents assum- 
ing values different from the mean field ones (note that 
the lower critical dimension of the Ising model is 1). 

Another interesting point which can be studied by con- 
sidering the Ising model on this network is the relation 
between the shortest path exponent 9 and the effective 
dimension of the network for 1 < S < 2. If one interprets 



9 as the inverse of some effective dimension, d e s, then, 
according to [15|. for some values of S > 1, d ff is actually 
greater than 4. In fact we have used the burning algo- 
rithm [24| to find out independently the shortest paths 
between all pairs of nodes in network of size N and found 
out that indeed 9 has values below 0.25 upto S ~ 1.2 (Fig. 
|5J . Our estimates of 9 are slighlty larger than those ob- 
tained in |l5j . which may be due to the smaller system 
size considered, but still we find 9 < 0.25 over a con- 
siderable range of values of S > 1. The question that 
naturally arises is whether the Ising model shows mean 
field behaviour here as the effective dimension is higher 
than the upper critical dimension, i.e., does it have a 
mean field behaviour even though the network is not a 
small world? 



III. ISING MODELS ON NETWORKS 

Ising models with long range neighbours have been 
considered in different contexts in many earlier studies. 
It is well known that the short ranged Ising model on a 
one dimensional lattice does not show any phase transi- 
tion while interaction with all the neighbours gives rise 
to a mean field like phase transition. The occurrence 
of the random lo ng rang ed bonds, also enables a phase 

transition HE EH Ea ISE EE 12a . 

In a analytical study of the Ising model on a WS type 
network, it was assumed that there are two nearest neigh- 
bours and additional long range neighbours occur ran- 
domly with probability p — p/N [l!j where p is a pa- 
rameter. Phase transition occurs here for p > 1 implying 
that the addition of one long range bond to each node, 
(on an average) could enhance a phase transition. In the 
numerical studies also, the addition of one random long 
range neighbour per site enables a phase transition |25j . 
The fact that the Ising model on these networks show 
a phase transition even when the additional long range 
connections are finite in number (i.e., O(N)) compared 
to that of the the infinite range model (where there are 
0(N 2 ) number of bonds) is an evidence of the key role 
played by the randomness in these networks (2(|. 

Another kind of long range Ising model has been con- 
sidered in which the interactions decrease with the dis- 
tance as a power law, i.e, the Hamiltonian is given by [28| 

H = —HijJijSiSj, (2) 

where ~ -jr ; r.y is the distance separating the spins at 

ij 

site i and j. In the fully connected model, phase transi- 
tions occur at the values p = 3/2 and p = 2; for p < 3/2, 
the system is meanfield like while for 3/2 < p < 2 it is like 
a finite dimensional lattice and for p > 2, it behaves like 
a one dimensional system with no phase transition. On 
a small world network, the interactions exist randomly 
among neighbours, maintaining the form of as above. 
However, there is no phase transition displayed by this 
model |2!j except at p = where it coincides with the 
WS model. 
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As mentioned earlier, the critical behaviour of Ising 
model on a WS addition type network is mean field type. 
The intriguing feature is the validity of the finite size scal- 
ing analysis in the mean field regime. Moreover, it has 
been observed |25j that the data collapse in the finite 
size scaling analysis can be achieved when \T — T c \, the 
deviation from the critical temperature is scaled by the 
factor N" 1 ^ in the scaling argument with D — 2. Inter- 
preting v as vd, where v is the critical exponent for the 
correlation function and d the effective dimensionality of 
the system, this could imply that one has effectively a 
system with d = 4 with the mean field value of v = 1/2. 
This effective dimension is identical to the upper critical 
dimension of the Ising model. 

Phase transition of the Isin g m odel has also been ob- 
served in scale-free networks |3(J, |31| where the transi- 
tion temperature diverges with the system size. Several 
other aspects of Ising models on networks have been con- 
sidered recently, for example the self-averaging prop erty 
of Ising model on net works |26j. quenching dynamics |32| 
etc., which are not directly related to the content of the 
present paper. 
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FIG. 3: Data collapse for different system sizes N = 
100, 150, 200, 250 for U, x, M and c, for 8 = 0.6 (p = vd). 



IV. ISING MODEL ON EUCLIDEAN 
NETWORK 

We have considered an Ising model on a finite chain of 
length N with periodic boundary condition. The nodes 
are assigned position coordinates 1, 2, . . . , N along the 
chain. Each node is connected to its two nearest neigh- 
bours. In order to generate the long range bonds, two 
nodes are selected randomly. If I > 1 is the distance sep- 
arating them, they are connected with probability P(l) 
as given by eq (1). The process is repeated N/2 times 
generating N long range bonds (each bond is counted 
twice so that there are N long range bonds in the sys- 
tem) . This ensures that there is one long range neighbour 
for each site on an average. As has been mentioned in the 
previous section, this is sufficient to achieve a phase tran- 
sition. For each realisation of the network, the dynamical 
evolution from a uniform state (all spins up) was allowed 
following a Metropolis algorithm for different tempera- 
tures T. 

We have computed the following quantities on this net- 
work after it reaches equilibrium: 

• Magnetisation per spin m = J^. Si/N 



Binder cumulant U = 1 — 



(r» 4 ) 
WWW' 



• Susceptibility per spin is calculated from the fluc- 
tuation of the order parameter: 



X 



N 
~K~bT 



[<m 2 } - (My 



where M is the total magnetisation, T is the tem- 
perature and Kb the Boltzmann's constant. 



• Specific heat per spin is calculated from the fluctu- 
ation of the energy of the system as 

where E is the energy of the system. 

A. Results 

From the intersection of U for different N in the plot 
of U vs T we estimate T c and a data collapse is obtained 
by plotting U vs (T - T^N 1 /". With this value of D, one 
can now estimate the exponents /3, 7 and a, the criti- 
cal exponents for the order parameter m, susceptibity x 
and specific heat c respectively using finite size scaling. 
Figures and 0] show the data collapse for these quanti- 
ties for S = 0.6 and S = 1.4. Similar results have been 
obtained for other values of 5. 



Transition temperature 

As expected we find a transition temperature T c which 
decreases with 8 (Fig. |SJ). T c varies very slowly in the 
region < 5 < 0.5 and much faster for higher values of 
S. 



Critical exponents 

From the finite size scaling analysis, we obtain the crit- 
ical exponents (3, 7 and a fFigs I(jl7fl . We find that v is 
equal to 2 for the entire range < S < 1. However, for 
5 > 1 we find that v decreases and then rises again. 
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FIG. 4: Data collapse for different system sizes N = 
50, 100, 150, 200 for U, x, M and c, for 8 = 1.4 (v = i^d). 
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FIG. 5: Variation of critical temperature T c with 5. Error 
in this and the next two figures are of the order of the size of 
the data points. 



The magnetisation exponent (3 is equal to 0.5 for 5 < 1, 
while for 8 > 1, it decreases with 5. The susceptibility 
exponent 7 is close to 1 for 6 < 1 and increases with 
6 beyond 5 = 1. The specific heat exponent a remains 
constant at a value close to 0.1 for 5 < 1 and appears 
to deviate from this value as 5 is made larger than 1. It 
becomes small as S approaches 2. 

The behaviour of all the quantities are consistent with 
the fact that for 5 < 1, there is a mean field behaviour 
of the Ising model (consistent with the available solution 
for the random long-range bond small world network |27j1 
and for 5 > 1, its behaviour is like a finite dimensional 
lattice with effective dimension greater than one. 

As has been argued for S = 0, here also one can inter- 
pret v as vd, where d is the effective dimension, such that 
it satisfies v — 2 — a. Since below d = 4 a is nonzero, v 
decreases from 2 at S > 1. Theoretically, the behaviour 
of a is non-monotonic as a function of dimensionality d 
as d varies from 4 to 1; a is zero for both d = 4 and d = 2 




FIG. 6: Variation of v and specific heat exponent a with 8. 




FIG. 7: Variation of magnetisation exponent /3 and suscep- 
tibility exponent 7 with 8. 



and this behaviour is reflected in the fact that v again 
rises to 2 at around 6 = 1.8. However, our result for a 
does not show this non-monotonic behaviour very accu- 
rately, the reason possibly being that the magnitude of a 
is small (O(0. 1)) even when it is non-zero and is difficult 
to estimate accurately in a numerical study. 

Mean field value for 7 is 1.0 which comes out to be 
slightly lower in our estimate and also the value of a is 
higher than the corresponding 0, but the scaling relation 
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a + 2/3 + 7 = 2 is obeyed with a high degree of accuracy. 
Although there are deviations of a and 7 from the mean 
field values, the fact remains that they show no variation 
over the region < S < 1, indicating that the critical 
behaviour for S < 1 is identical to that for 6 = 0, which 
is known to be mean field like. 



B. Effective dimension 

Our results confirm that the network behaves as a finite 
dimensional system in the region 1 < 6 < 2. However, 
the finite size scaling analysis of various quantities in the 
Ising model does not allow the estimation of d and v inde- 
pendently. It may be expected that the behaviour of the 
average shortest path should give us some idea about the 
effective dimension. However, as we argue below, it is not 
possible. The estimate of the exact average shortest path 
is consistent with the mean field picture for 5 < 1 as it 
shows a logarithmic increase with N (small world effect). 
As mentioned before, for 6 > 1, it has a power law in- 
crease N e and we have a region where 9 < 0.25. So if one 
interpretes 9 as a inverse dimension, the dimensionality is 
greater than 4 here. Even then, the Ising model does not 
show mean field behaviour here. From this we conclude 
that although there is an effective dimension coming out 
from the topological behaviour of the network, it does 
not exactly act as the spatial dimension as far as physics 
of Ising model is concerned. 

TABLE I: Table showing the effective dimension and compar- 
ative values of vd from e expansion (ee) and v from finite size 
scaling (fss) analyses. 



S V e g 


deS 


vd (ee) 


v (fss) 


1.2 0.56 


3.34 


1.85 


1.85 


1.4 0.59 


2.90 


1.72 


1.75 


1.6 0.61 


2.74 


1.66 


1.9 


1.8 0.72 


1.36 


0.98 


2.0 



We adopt a different approach to estimate v and d 
which works well for S close to 1. The behaviour of the 
critical exponents is fully consistent with the fact that the 
effective dimensionality d decreases as S is made larger 
than 1 (e.g., the value of 7 increases while that of [3 
decreases). In order to find out d, we use the epsilon 
expansions of a, and 7 to first order, each of which 
gives an independent estimate of e = 4 — d: 

e = 3(1- 2/3) 
e = 6(7-l) 
e = 6a. 

Ideally, all three estimates should give similar value of e, 
but here we have used the expansions upto linear terms 
only and the estimates may not agree very well (espe- 
cially for the estimate from a and to some extent, that 
from 7). It is not very convenient to use the expansions 



to higher degrees also. To take care of this fact, we take 
(e), the average of the above three estimates of e, and use 
it to find out the value of v from the equation 

Thus an estimate of vd is obtained which can be com- 
pared to v obtained from the finite size scaling. We see 
that when S — 1 is small, the agreement is nice (Table 1). 
Naturally, when 5 is considerably away from 1, the effec- 
tive dimensionality decreases and the epsilon expansions 
are not much useful here, particularly for a and 7. This 
is reflected in the mismatch of the estimates from epsilon 
expansion and the finite size scaling analysis for S = 1.6 
and S = 1.8. 



V. SUMMARY AND CONCLUSIONS 

We have investigated the behaviour of the Ising model 
on a one dimensional network in which there are links 
to random long range neighbours existing with a proba- 
bility which varies inversely with the distance along the 
original chain following eq The effective dimension- 
ality of the network deviates from 1 continuously as S is 
made smaller than 2 and a finite temperature phase tran- 
sition is observed. The transition temperature increases 
as 5 is made smaller and the critical exponents vary with 
5 as the dimensionality changes. The results show that 
there is a mean field behaviour for S < 1 and a finite 
dimensional behaviour for 1 < S < 2. In the mean field 
regime, the finite size scaling analysis works with an ef- 
fective dimensionality four, the upper critical dimension 
of the Ising model. The present study confirms the con- 
clusion made in that the small world behaviour of 
the underlaying network exists for 5 < 1 only. 

Our results also show that the effective dimensionality 
cannot be simply extracted from the behaviour of the 
shortest paths on this network although the latter shows 
a power law behaviour with the number of nodes. We 
have used an alternative method to estimate the effective 
dimensionality which works quite well for 5 close to 1. 

The present study was conducted to address certain 
issues related to the topology of a network but the study 
of phase transition of Ising model on networks has its 
own inherent interest as the Ising model has got wide 
applications in many interdisciplinary areas. For exam- 
ple, the study and nature of phases in such models gives 
us important insight into social opinion dynamics [33J in 
closely knit populations, where the networks look very 
similar to the models discussed here, and opinions being 
modelled as states of spin vectors at each site. 
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